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Q ' Abstract 
(N ■ 

, We consider a totally asymmetric exclusion process on the positive half-line. 

\ When particles enter in the system according to a Poisson source, Liggett has 

computed all the limit distributions when the initial distribution has an asymptotic 
density. In this paper we consider systems for which particles enter at the boundary 
according to a complex mechanism depending on the current configuration in a finite 
neighborhood of the origin. For this kind of models, we prove a strong law of large 
numbers for the number of particles entered in the system at a given time. Our 
main tool is a new representation of the model as a multi-type particle system with 
' infinitely many particle types. 



1 Introduction 

^ \ The simple exclusion process r], = {r]t)t>o on a countable space S, with random walk 
\^ \ kernel is a continuous time Markov process on X := {0, 1}'^. For a configuration 
■ G X, we say that the site x is occupied (by a particle) if ri{x) = 1, and is empty if 

r]{x) = 0. A particle "tries" to move from an occupied site x to an empty site y at rate 
Q . p{x,y), or in an equivalent way, waits for an exponential time of parameter 1 and then 
0\ ! chooses randomly a site y with probability p{x,y) and "tries" to jump on y. When the 
^ ' site y is already occupied, the jump is canceled and the particle stays at x. In this way, 
there is always at most one particle at a given site. Formally, the exclusion process 77. is 
^ ' defined as the Feller process with generator 



(if('l) E P(^, VH^) (1 - >)(»)) Ifil,.,) - fil)] , (1) 

x,y£S 

for all cylinder functions /, where 

{r]{y) if z = x, 
ri{x) ifz = y, 
ri{z) otherwise. 

A natural question is to describe the set of invariant probability measures 1 which is the 
set of probability measures ^ on S such that, if r/o ~ then for alH > 0, r/j ~ /i. These 
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measures are characterized by the equations: 

j nffiidv) = 0, 

for any cyhnder functions / (see e.g. [5| for a review). In the case S = the set of 
translation invariant stationary measures is exactly the set of convex combinations of 
product Bernoulli measures on Z (see 

In this paper, we consider the case S := Z"^, or S := {!,...,/?} with R > 2, and 
p{x, x + 1) := 1, i.e., the totally asymmetric nearest neighbor case. In Z^, one has to add 
some boundary mechanism to make the model non trivial. The simplest way to do this 
is to add a particle reservoir at site with a certain density A > 0. This means that a 
new particle is created at site 1 according to a Poisson process with rate A when this site 
is empty. We call the model on the TASEP(A), and we denote by Qx its generator 
and by S\{t) its semi-group: 

oo 

njir^) = A (1 - r^(l)) [/(r^i) - /(r^)] + J] r^{x) (1 - r/(a; + 1)) [f{v.,.+i) - /(r/)] , (2) 

x=l 

for all cylinder functions /, where 

.,(.):=( ^f^=^' 
1^ ri{z) otherwise. 

In (I2]) we see two parts for the generator: one is due to the boundary mechanism and 
we will call it the boundary part; the other one, which has the form given by ([T]) for 
S = Z^, is due to the exclusion process and we will call it the bulk part. In the same way, 
if the state space is {1, ... , R}, then we add a particle reservoir with density p G [0, 1] at 
site R+ I, which means that when the site R is occupied, then the particle in this site 
disappears with rate 1 — p. 

Let us introduce some notation. In the following, we denote by z/'^ the product measure 
on Z^ with density A and by r the shift, r acts on configurations ?7 G X by 

Tri{x) := ri{x + 1), Vx G Z^, 
on functions / : X — > M by 

rfiv) ■■=fiTr]),yrjeX, 
and on measures /i on X by 

j fdTp:= J Tfdp,"^/ e L\fx). 
For a measure /x on S* and / G L^{p), we will denote (/)^ = / fdp. 

We are interested in the limit behavior of the distribution at time t. For this model, 
we have a good understanding about what happens at equilibrium. Indeed, Liggett has 
shown in [3] the following ergodic theorem, which gives the limit measure for an initial 
measure with a product form and an asymptotic density: 
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Theorem 1.1 (Liggett [3]). Let n be a product measure on for which 
p := lim^^oo {v{x))n exists. 

1 I A^n: if P > h (bulk dominated), 
IfX>- then lim 7iSx(t) = { / ' ]\ ^' 

2 t^oo 1^/^15 ¥ P ^ 2 i^o-xiTnum current). 

1 U^, if p>l-\ (bulk dominated), 
If X < - then lim nSxlt) = s ^ 

2 t^oo [i^ , if p ^ I — ^ (boundary dominated) , 

where the p^ 's, for p > |, are stationary measures and asymptotically product with density 
p, i.e., lim^^ac 'T^ f'p = ^'^ (in a weak sense with test functions f G C{X,'R)). We also 
have pI = u^. 



To describe the set of invariant probabihty measures in the cases S = Z and 5* = , 
he uses extensively that the product Bernoulh measures are invariant and for these 
mesures one can make exphcit computations. To break the product form of the invariant 
measures, one can consider a boundary mechanism which is not Poisson. We hmit our- 
selves to finite range boundary mechanism, i.e., for which there exist some R G such 
that the boundary part of the generator vanishes on every cylinder function with support 
in {R + 1, . . .}. This idea was first introduced by Grosskinsky in chapter 3 of [1] where 
he defined the following Feller process: 

^fiv) ■■= E ^(^) (1 - + 1)) [f(v.,.+i) - fiv)] 

(3) 

+ E k.i'\s,=iW^^\^SR)-f{r^)\. 

for all cylinder functions / where S^^ := {1, . . . , i?}, := {0, 1} ^ and ^' U ti^Sr is the 
natural concatenation of configurations on and on ^S^. In ((H]), the first sum describes 
the bulk part and the second sum the boundary part of the dynamic. 

The reason for which we only treat the finite range case is that when we are not in 
this case, pathological things can occur. For example, consider the following dynamic 
with non local boundary mechanism: Define the asymptotic density of a configuration 
G X by p(r/) := liminf^^oo \ ^YTi=\ v{i)'^ 'we consider now a TASEP on Z^ for which the 
rate of apparition of a particle in site 1 is p{ri) where t] is the current configuration. More 
formally, the boundary part of the generator is 

p(r/)(l-r^(l))[/(r/i)-/(r/)]. 

In this example, every mixture of product Bernoulli measures is invariant for the process. 

For this generalized boundary mechanism, we will not have an exact solution as for 
the TASEP (A). Our approach is to study the number of particles entered in the system 
up to time t. We will see that it grows linearly with an almost sure speed equal to the 
stationary current joo- Define poo as the root of p(l — p) = joo in [0, l/2[. We think that 
the process has a stationary measure which is asymptotically product with density poo 
but we are still unable to prove it. 
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The rest of the paper is organized as follow: in section [2] we give a construction of 
the process defined above using a graphical representation similar to that introduced by 
Harris [2]. We also introduce the basic coupling technique which is the main tool used 
in the paper; in section Owe give some general results on the asymptotic behavior of the 
TASEP with complex boundary mechanism; finally, in section H] we study a particular 
example: take a TASEP (A) on and add a source (independent of everything) with 
density e > which is activated only when site 2 is occupied. For this model, let Nt be the 
number of particles in the system at time t when we start from the empty configuration. 
Then the main result of this paper is the following strong law of large numbers: 

Theorem 1.2. Almost surely, 



where p{X) is a positive constant (depending only on X) for which we give a natural 
probabilistic interpretation. 

It should be noted that this particular choice of boundary mechanism is rather ar- 
bitrary, and that our method is robust enough to be used in a much larger generality. 
However, the notations which would be needed would be much more tedious, while pro- 
viding very little additional insight into the model — so we choose to limit ourselves to 
one representative case. 
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2 The Harris construction 

We will use the method developed by Harris [2] to construct our process. Let 



be a family of independent Poisson point processes on constructed on the same prob- 
ability space (r,jF, P), such that the rate of processes indexed by Z+ is 1 and the rate 
of the process indexed by (r/, is > 0. By discarding a P-nuU set, we may assume 
that 

each poisson point process in J\f has only finitely many jump times in every , , 
bounded interval [0, T], and no two distinct processes have a jump in common. 



t— >oo t 



lim — = A(l - A) + A(l - A)p(A)e + o(e) 




We denote 



U 
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Fix T > and r] e X. The process {ilt)o<t<T starting from tj is now constructed as 
follows. Consider the following subgraph of 

Gt := + : X > i?,A4n [0,T] ^ 0} a; + 1} : x G {0, . . . , i? - 

It is easy to see that every connected component of Gt is almost surely finite. Let Fq be 
the subset of F such that (jlj) and the above condition hold for all T > 0. Then we have 
Fq G and P [Fq] = 1. We consider now only u G Fq. For every connected component 
C of Gt, the set (U^^gc-Mc) f] [0'^] finite so its elements can be ordered chronologically 
Ti < . . . < Tn and we need only to describe the action of each of them. We start with the 
configuration r]: 

r]t{x) := i]{x) 

for all X G C and < t < ri. 

Suppose that the process is constructed on C for < t < and A; G {1, . . . , n}. Then: 

• if Tk G A/'^,^' and if ^r^r |5jj = ^ then r]r^\Sji ^' Vrk {^) '■= Vr- (^) ^ ^ C\Sr, 

• if Tk G A/'^.g' and if Vt^\Sr ^ then rir^{x) := rj^-{x) for all a; G C, 

• if Tfc G Ux and r] - {x) ( 1 - r] - (x + 1) ) = 1 then 77^-. := ( - ) on C, 

• if Tfc G A/"!: and 77^- (x) ^1 — rj^- (x + 1)^ 7^ 1 then 77^-^. := 77^- on C, 

Finally, we put rjt := 7/^^. on C for t^ < t < t^+i ii k < n and for Tn < t < T ii k = n. We 
make the same construction on every connected component of Gt and then let T go to 
infinity to get the process {r]t)t>o ^'^^ every lu G Fq. 

The usefulness of such a construction is that, using the same Harris process, we can 
construct two or more realizations of the process on the same probability space starting 
from different initial configurations. We will refer to this coupling as the basic coupling. 

3 The attractive case 

In this section, we consider the process with the generator: 

(5) 

+ E 4,e'W«[/(e'U7/|.5j-/(r7)], 
and we assume the process is attractive. 
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3.1 The stationary measure 

Proposition 3.1. Suppose that the process is monotone (or attractive). Then, starting 
from the empty configuration, the measure at time t of the process, say fit, is stochastically 
increasing and converges to a measure fi^o G X which is the smallest invariant measure 
of the process. Furthermore, fi^o G Te a'l^d /loo is ergodic. 

Proof. Let < s < t. We have 5q -< fit-s- Thus by attractivity of the process, we have 
6oS{s) -< fit-sS{s), i.e., -< fif. Hence, by monotonicity, fit converges weakly to an 
invariant measure /loo- 

For all 1/ G X, we have 6o -< u, which implies that fit -< v for all t > 0, and then 
fioo -< Assume now that fi^o = Az/i + (1 — A) z/2, with z/i, z/2 G X and A G ]0, 1[. We 
have fioo = Az/i + (1 — A) z/2 >- /^oo, thus 1/1 = 1/2 = fioo and fioo is extremal. Finally, by 
Theorem B52 of [5], fioo is also ergodic. □ 

Proposition 3.2. fioo is stochastically dominated by the product Bernoulli measure 
with density j. 

Proof Define TV' := {M^,xe Z+), where := K+r- Then Af' defines a TASEP (^t) 
on with rate 1 of particle apparition in 1. By theorem 11.11 starting from the empty 
configuration, the distribution at time t converges to z/5. In this coupling, we have 
^t{x) > rit{x + R) almost surely for all t > and x > 1. Thus the restriction of fioo to 
{R + 1, -R + 2, . . .} is stochasticaly dominated by 1/2. □ 

3.2 Asymptotic measures 

Let us extend the measure fioo on a measure on {0, 1}^ by 
7^oo(^) ■■= fiooiv ^ X -.f] e A}, 

w s if X > 1, 

where rj[x) := < 

I otherwise, 

for all A in the cylindric field of {0, 1}^. We will still denote this measure fioo- Let 
flk '■= T^f'oo and consider any weak limit fioo of this sequence: 

lim flk, = fioo- 

i—too 

Proposition 3.3. The measure fioo is a translation invariant stationary measure for the 
TASEP on Z. Consequently, it is a mixture of product Bernoulli measures, i.e., there 
exists a probability measure a on [0, 1] such that 

Jo 

Proof. Let fi^ be the generator of the TASEP on Z. For any cylinder function / : 
{0, 1}^ — > M, let X G Z^ large enough such that supp r^f C {R + 1, R + 2, . . .}, where 
suppf is the support of /. Thus r^f could be considerate has a function on Zl and we 
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can apply the generator to this function. We get fir^/ = ^I'^r^f for all y > x. But it 
is easy to see that fi^ and r commute, thus we have 

= J n^fUdv)- 

Hence for i large enough, = 0, which implies that (fi*^/)-^ = 0. This is true for 

arbitrary / thus /too is invariant for the TASEP on Z. We know that for this model we 
have Xg = A e [0, 1]} U {i^^, n G Z}, where Vn = 't"'^o aiid is the Dirac measure of 
the configuration for which all the sites a; > are occupied and all the sites x < are 
empty (see [1]). Using Proposition 13. 2^ /ioo is stochastically dominated by thus fioo is 
translation invariant and is a mixture of product Bernoulli measures. □ 



3.3 A strong law of large numbers 

Let n be an invariant and ergodic measure for the process with the generator given by 
([5]). Fix ^0, ^ and ^' three configurations on Sr and consider 

Nit) ■.= u^fu'r^It), 

where It '■= {s E [0, t] : r]s\Sii = ^o}- We will show a strong law of large numbers for N(t) 
which will be useful in the sequel. 

Proposition 3.4. // the process starts from fi and if ^' ^ C,o, then almost surely: 

lim ^ = [ri E X : risisj, = ^o} . 

Proof. Let 

*J 

and 

i){t) := inf {s > : T, = t}. 

Since /i is ergodic, — ^ {rj E X : r]s\s„ = almost surely. Let I := {t > : msj^ = 

^o}- ^ : IR+ — > / is a one to one map so we can define ^A := ip^^{Af^^^' fl /) and N'{t) 
the associated counting process. We have now that N'{t) = N{tlj{t)) almost surely and 
is a Poisson point process with parameter ci^,^'. Indeed, for all t > 0, the next point 
in Ai after t is equal to t plus an exponential variable with parameter d^^^i, and after 
this point, the process is independent of the past (see Figure [T]). Then, it is easy to 
see that for every e > 0, ipiTt) < t < ip{Tt + e). Since N{t) is non-decreasing, we get 
N'{Tt) < N{t) < N'iTt + e). Consequently: 

N'{Tt) Tt ^Nitl^ N'jTt + e)Tt + e 
Tt t - t - Tt + e t ' 

Since both sides converge to d^^^in[rj E X : ?7s|s^ = ^o} almost surely, it leads to the 
conclusion. □ 
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7N 7N 7k TK 7i\ 7M\ 7K 7K TK 7K 

I r I B ^ B 

,,.'(*) 

Figure 1: On the time interval [0,ip{t)] we see the set I^(^t) in gi"^y- 
The total length of the grey part is t. The stars are points of the 
process A/'^,^'- In this example, N'{t) = 5. 



4 A particular case and the Multi-Species model 

In this section, we are interested in a particular case of TASEP with a complex boundary 
mechanism: let A, e > such that A + e < |. Particles are created at site 1 with rate 
A + e7]{2), where rj is the current configuration and the bulk dynamic is the one of the 
TASEP. This model has a generator given by: 

^fiv) = Yl ^(^) (1 - + 1)) if(v.,.+i) - fiv)] 

x&x (6) 
x(l-r/(l))(A + 6r/(2))[/(r/i)-/(^)], 

for all cylinder functions / on X. As it is explained in the introduction, the choice of the 
model is rather arbitrary, and the methods that we use are quite robust (at least as long 
as the system can be dominated by a product Bernoulli measure of intensity lower than 
1/2 — which is indeed the case here). 



A 




Figure 2: Particles enter with additional rate e when the site 2 is 
occupied. 

In this model, the range of the boundary mechanism is i? = 2. The hypothesis e > 
implies that the process is monotone, thus we can define the smallest stationary measure 
/^oo = Z^oo (A, e) of the model. Using the Harris representation, we can couple this process 
with 7]^, a TASEP(A), and r/^+^ a TASEP(A + e), such that if r]o < Vo < Vo^' then for 
alH > 0, rj^ < rjt < ^t^^ ■ This proves that -< fi^o -< ^^^^ and then -< floo -< z^^"^^- 



4.1 Some estimates about the flux 

Here we will see another way to see the process with the generator given by (Q. For any 
i > 1, let 

Xi := {oo, 1, . . .,i}^+ . 
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We define 



+ l]el^(i)>j+il^(2)=j-i [/ iVj^i) - fiv)] 

(7) 

+ el,,(l)=oolr,(2)=i [/ in^^l) - fiv)] 



x=l 

for all cylinder function / : M, where 



j if X = 1, 
r]{x) otherwise, 



for j G {1, . . 
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Figure 3: First class particles enter with rate A whatever is the 
configuration in {2, 3, . . .} and second class particles enter with rate 
e if the site 2 is occupied by a first class particle. Particles in black 
are indistinguishable particles (their class has no influence on the 
rate of the source in the current configuration) . 

We fix 2 > 2 for the sequel. The new description is described in Figure [S] and in the 
following. We put the particles into a certain number of classes. For a configuration 
rj G Xi and for a site x G Z^, the number rj{x) designates the class of the particle at 
site X if it exists, i.e., if ri{x) ^ oo, and is equal to oo if the site is empty. We use here 
another notation for empty sites because it allows us to have a simpler expression for the 
generator and we can also interpret holes as particles of class infinity. The evolution is 
the same as before, except that if a particle of the fc-th class (or of type /c) attempts to 
jump on a site occupied by a particle of the j-th class (or of type j), then it is not allowed 
to do so if > j, and the particles exchange positions if < j. We say that a particle of 
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class A; G {1,2,.. .} has priority on all particles of classes greater than k. In this way, a 
particle of type k behaves as a hole for particles of type j < k. 

Now we will explain how we affect classes to the particles. First class particles enter 
in the system (in the site 1) at rate A. As they have priority on other particles, they 
are not affected by them, so the process of first class particles is simply a TASEP(A) on 
Z^. Next, particles of class 2 < j < i — 1 enter in the system with rate e, if the site 2 
is occupied by a particle of class j — 1 and with rate otherwise. Finally, particles of 
class i enter in the system with rate e if the site 2 is occupied by a particle of class i — 1 
or i and with rate otherwise. For each configuration of the system, at most 2 types of 
particles are allowed to enter in the system. We can also remark that if we consider the 
process consisting with particles of class 1, . . . ,i, then it has the generator given by 

In terms of Harris system, we define J\f the set of the following independent Poisson 
point processes on M.^: let (A^, a: > 1) be Poisson point processes of rate 1; let (Ar/,j > l) 
be Poisson point processes of rate A for AAf and of rate e for the others. The mechanism 
is then the following: if t > is a jumping time of A/^. and if at time t~ we have 7]{x) ^ oo 
and ri{x + 1) > ri{x) (i.e., there is a particle at x and it has higher priority than the one 
at X + 1 if it exists), then the particle at x jumps to x + 1 (and the particle at x + 1 jumps 
to X if it exists); if t > is a jumping time of N'l and if at time t~ we have ri{l) > 2, then 
a first class particle appears at site 1; if t > is a jumping time of Afj with 2 < j < i — 1 
and if at time t~ we have 77(1) > j + 1 and 77(2) = j — 1, then a j-particle appears at 
site 1; finally, if t > is a jumping time of Afji' and if at time t~ we have ri{l) = 00 and 
77(2) E {i — I, i}, then a i-particle appears at site 1. 

We denote by S^'^{t) the semi-group corresponding to the generator fi^*) and by 
iVt''^)t>o the process of the j-th class particles for j = 1, . . . ,i. It is easy to see that 
the distribution of the process {^'t^)t>o for j G {1, . . . , i — 1} does not depend on i. The 
process is attractive, thus we can define /x^^ as the weak limit of 5oo>S''^*^(t). As in Propo- 
sition 13.11 this measure is extremal, ergodic and the smallest invariant measure of the 
system. For all 1 < j < we note ffj.^^ := 'Yl^k=i'nt^'^ ■ Remark that the process {f]f')t>o 
is exactly the process that we want to study, i.e., it has the generator given by Fur- 
thermore, for all j G {1, . . . , the distribution of the process {f][^^)t>o is the same for all 
i. 

In order to compare the processes (f/f ^^)t>o and (r7f'')t>o, we need to control the 
number of particles of a given type in the system at a given time. Let N^^^ be the number 
of j-particles entered in the system between times and t, and define 



for all j G {1, . . . , i}. 

T^^^ is the t ime spent by j-particles in site 2 during [0,t], and Tt is the length of 
the subset of [0, t] for which 2-particles can enter in 1 with rate e (excepted if the site 1 




and 
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is already occupied by another 2-particle). The foUowing lemma says that we have an 
uniform control on the total time spent by a particular particle of type > 2 at site 2. Let 
Tfcj be the total time spent in site 2 by the k-th particle of type j > 2 entered in the 
system. 

Lemma 4.1. There exists a constant Ca G ]0, +oo[ such that for all k > 1 and all j >2 

we have 

E [n,] < c 

Proof. Let Ef be the event that, between times t and t + 1, a first class particle enters 
(or tries to enter) in the system, then it jumps, if it is possible, to the site 2 and finally 
another first class particle enters in the system. We also assume that in Et there are no 
other jumping times for A/i, A2 and between and t. In particular, if Et occurs and 
if there was a particle of type > 2 in site 2 at time t, then it has disappeared at time 
t + 1. We have that P [Ef] = g(A) does not depend on t and g(A) > 0. 

On the event {T^j > t}, there exists a time r such that the A;-th particle of type j 
is at the site 2 and it has spent exactly a time t in this site between and r. We have 
Er C {Tkj <t+l}. Hence 

P [Erinj >t]<P [Tk,j < t + l|Tfc,,- > t] . (8) 

But r is a stopping time for the Markov process {r]f \ I = 1, . . . ,j)t>o and the event Er 
depends only on the poisson processes of the Harris system for times between r and r + l, 
so, conditionally at {r < 00}, E^- has the same law as Eq by the strong Markov property. 
Hence the left hand side of ((H]) is equal to q{X). Finally, we have 

P [Tk,j > t + 1] < (1 - g(A)) P [Tk,j > t] . 

The last inequality implies that there exist some deterministic positive constants ai,a2, 
depending only on A, such that almost surely and for all t > we have 

P [Tk,j >t]< a^e-'^'K 

The result follows with Cx := Jq°° aie~°''^^dt. □ 
Finally, the following theorem gives the estimates that we need: 

rpU) 

Theorem 4.2. For all 1 < j < i, converges almost surely to a deterministic value if 

(k) 1 

the process starts under fioo for all k > i. Furthermore, for all e < where C\ is as 
in Lemma \4-1\ we have 

limsup — - — < Cj-ie-' ^, lim — — < cje^ ^, 
for 1 < j < — 1, and 

limsup < 2ci-ie''\ lim — < 2cie'-\ 

t^oo t t^oo t 



where {cj)j=i^,„^i are constants (depending only on X) such that cq := A (1 — A) and Cj 

X 
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Proof. We have seen that every /i^-* is stationary and ergodic, so by the ergodic theorem, 
we have almost surely 

^ f^^^Hv ^ X, : r^{2) = j} (9) 



t t- 



>oo 



and 



T {veXk-. Vil) > 2, r/(2) = 1} . (10) 

Since the distribution of the first class particles is u under every /x^o , the right hand 
side of ([9]) is A if j = 1 and the right hand side of ( fTOj) is A (1 — A). Using Proposition 
13.41 N^^^ /t converges to A (1 — A) almost surely. 

Let 

Mf ■.= ^{se M', n [0, t] : r^«(2) (l - r^«(l)) = l} . 
Then almost surely A^^^ < and applying Proposition 13.41 

lim sup ^— < eA (1 - A) = lim — ^. (11) 

Now, we need to find an upper bound for limt^oo • First, we can remark that 
can be decomposed into two parts: the time spent by initial second class particles T^^ 
plus the time spent by the new second class particles T^g'* in site 2. But, since t/^"* is 

bounded by a random variable almost surely finite, it is sufficient to study limt^oo 
As we have seen previously, jji'S -< v^^'^. The idea is that since we know the number of 
second class particles created up to time t, it is sufficient to bound the time spent in site 2 
by one of them in the environment v^^'^ where it is slower. But there are some difficulties. 
For example, at the moment where a second class particle is created, the environment in 
{2, 3, . . .} is not dominated anymore by a Bernoulli product with density A + e because 
we know that a first class particle has to be in site 2. To avoid this problem, we will use 
the following fact: if a particle of a class different than 1 is at the site 2 at time t then it 
has a positive probability (depending only on A) to be out of the system at time t + 1. 
This implies Lemma [4.11 which says: 

E [Ti] < Cx, (12) 

where Ti is the total time spent by the l-th second class particle at site 2 and Cx is a 
constant. Take any /? > eA (1 - A) and r := inf |t > : Vs > t, nP < /5s|. We have 
that T is almost surely finite by (ITT]) and 

-^l{r<t} < -^Ta|^<j} < -^ril|^<4}. (13) 

i=l i=l 
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Taking expectation in both sides, it leads to 



E 



(2) 



t,2 



t 



l{r<t} 



I3t 



i=l 



Hence, by dominated convergence we have almost surely 



(2) 



(2) 



lim — — = lim — ^ = lim E 

t— >oo t t^oo t t—*oo 



(2) 



t,2 



L{r<t} 



</3Cx. 



The above inequality is true for all /3 > eA (1 — A), thus we also have 



(14) 



(15) 



lim 



(2) 



< eA (1 - A) C 



t—too t 

Let now C2 := CxCi and by induction, using exactly the same arguments, we have for all 
1 < j < ^ - 1: 

(i) 



lim sup — - — < Cj-ie-' ^, 

t— >oo t 



and 



lim 

t—*oo t 



where Cj := C{ A (1 — A). 

Finally, let a := lim sup j_ 
ffT^ . we get: 

lim — <aCx. 

t— >oo t 



► oo t 



Doing the same computation as in (fT3|) . ffT^ and 



Consequently, 



rp{i-^) _|_ rp{i) 
t^OO t 



lim 



which implies as in f llip : 



hm sup 

t— >oo 



a < (ci-ie' ^ + aCx)e. 



Since e < we have a < 2ci_ie* and 



lim 



(0 



□ 



Now, let A^j* and A^^*'''' be the number of particles entered in the system between 
and t for the processes {fit^~^^)t>o and {f]f^)t>o- We deduce from the above theorem, by 
summing estimates, that limsupj^^^ ' ~^ * = O (e*~^). 



r(0 
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4.2 The asymptotic flux at the flrst order 

In this section, we consider the particle system with the generator given by (JTj) for i = 3. 
In the previous section we have seen that if we want to compute the limit, at the order e, 
of the total number of new particles at time t divided by t, then it is enough to compute 
this limit only for the first and second class particles. In other words, if nJ''"^^ denotes 
the number of j-particles entered between and t which are still alive at time t, then: 

^T,(l)^ T,(2)^ T,(3) iVf'«+iVp) 

limsup = limsup h o(e), 

t—^oo t t— >oo t 

= A(l - A) + limsup — h o(e). 

T ("2) 

We will use the notation A^^ rather than A'"^ for the number of second class particles 
because there will be no possible confusion. The aim of this section is to prove a law 
of large numbers for A^^ and to compute the limit at the order e. Let c > such that 
A + c < I and e G [0,c]. Let us introduce some notations. Let Xj(t) be the position 
at time t of the i-th second class particle entered in the system, with the convention 
Xi{t) := if the corresponding particle either is not yet born or is died at time t. We 
define 

:= inf {t > : Xi{t) = 1} , 7^^ := sup {t > : Xi(t) > 1} , 

Si{t) := lxi(t)>i and Si := li;s=oo- 



# • • o 

I \ \ \ \ \ \ \ 

• o • 

I \ \ \ \ \ \ \ 

Figure 4: First class particles, in black, enter with rate A whatever 
is the configuration in {2, 3, . . .} and second class particles, in grey, 
enter with rate e if the site 2 is occupied by a first class particle. 

At time we start with the invariant measure fi^^ and we have -< fi^^ ~< u^^^. 
Since the dynamic is monotone, we can couple our process r/. with a TASEP(A) r^*"-^ and 
a TASEP(A + c) 77*"^ such that for all t > 0, we have rjl"'^ < Tjt < almost surely. We 
define a new particle system with state space {0, 1, (2,-i)j>i} and generator: 

^yfiv) ■= ^'^^{im ifiVi-^i) - fiv)) + el,,(i)^i,,,(2)=i {f{V2^i) - fiv)) , 

00 

+ lT7(a;)^0,r,(x+l)7^1 {f{Vx,x+l) " f (v)) , 

x=l 

for all cylinder function /, where 

\ r](z) otherwise. 
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and 




if z = 1 and ri{l) = 0, 
r]2^,{z):={ (2,2 + 1) ifz = landr7(l) = (2,z), 

otherwise. 

This particle system has the following description: particles of type 1 are first class 
particles and they have the generator Q^; particles of type 2 are second class particles, 
they enter with rate e if there is a first class particle at site 2 and they do not interact with 
any other second class particle. To link our process to the above process, we proceed as 
follows: at each time for which a second class particle enters in the real system, we add a 
supplementary particle to each of the processes r/*"-^ and r]^^^ if it is possible and we denote 
by Xl^-^{t) and X-^^{t) the corresponding trajectories of the particles associated to Xj. 
What we mean by supplementary particle is that these particles does not belong to the 
processes 77*"-^ and t^'^'^p, i.e., they move as second class particles and they do not interact 
with each other. In particular, a given site in one of these two processes can either 
be empty, or contains one first class particle, or contains one or more supplementary 
particles. Then the process t]'^"'-^ plus supplementary particles has the generator Clx and 
the process 77''"^ plus supplementary particles has the generator We have almost 

surely X-^^{t) < Xi{t) < X^^^ {t) for all t > 0. We define analogously the quantities 
^m/^ iV;"P, 7;^'-/, T^^'^-P, etc. 

Consider the following initial configuration: at time we put on {1, 2} first class 
particles according to (the product Bernoulli measure with density A), and we put 
one first class particle in site 2 and one second class particle in site 1. This is exactly 
the distribution of the configuration rffl for all z > 1. Then first class particles enter 
in site 1 with rate A and they have priority on the second class particle. Two cases can 
occur: either the second class particle survives, or it dies. Let p (A) be the probability 
that the second class particle survives, p is a non- increasing function, p(0) = 1, p(|) = 
and p (A) > for all A < |. Indeed, for the last point, it can be shown that if the second 
class particle survives, then it has a positive speed 1 — 2A (see e.g. [7]). We have by 



= p{\) and P [Sf"^ = 1] = p('^ + c) for alH > 1. Consequently 
<p{\ + c). 



construction P S"*"-^ = 1 

we have p{\) <Y* [Si = 1 

In this section we prove the following law of large numbers: 

1 Nt 

Theorem 4.3. Almost surely, lim - lim — = A (1 — A)p(A). 

eiO e t^oo t 

With the discussion at the beginning of the section. Theorem 11.21 follows. 

The idea is the following: when e is very small, second class particles do not interact 
before they are very far from the left boundary and if a second class particle is far enough 
from this boundary, then it survives with high probability. In other words, the effect on 
Nt of interaction goes to with e. The first part of the proof is to find estimates for the 
process without interaction and to prove the theorem in this case. Next, we will show, 
for the "true" process, that if two second class particles meet, then they both survive 
with a probability going to 1 as e goes to 0; this implies the theorem. 
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4.2.1 The process without interaction 

Consider a family {J^x) q^^^^i of Poisson point processes such that the parameter of Af^ 
is A and for all < A < yU < i, A/"^ C A/"^ and A/'^\A/'| is independent of X^- Take also a 
family (?7o)o<A<i i^iitial configurations such that 77o (2) (l — Voi^)) = 1 ^ill A e [O, | [, 
the distribution of rj^ on {3,4, . . .} is u'^, and for all x > 3 and all < A < /i < |, 
'7o(^) ^ Voi^) almost surely. Then using the same Poisson point processes {J\fx,x > 1) 
for the bulk dynamic we construct, as in section [21 the family of TASEP {i]^)q<x^i. such 
that 1]^ is a TASEP (A) and for alH > and all < A < < i, r/^^ < r]'^ almost surely. 
At time we add a second class particle in site 1 to each of these processes and we 
denote Xx{t) the position at time t of the particle in the process 77'*' (with the convention 
Xx{t) = if the particle is died). We define 

cA ._ -1 

• ^Xx survives; 

and 

:= M {t > : Xx{t) = x} , 
for all X > 1. 

We start with an intuitive lemma which will be useful to propagate results from the 
process without interaction to the true process. 

Lemma 4.4. The function p : [0, 1] [0, 1] is continuous. 

Proof. We start by proving the right continuity of p. Since p{X) = for A > |, it is 
sufficient to prove it on [O, | [ . Let < A < |, e' > and < c < A — |. There exists 
some X > 1 such that 

P [S^+^ = 0\T.^+^ < 00] < e'. 

Indeed, if M := max {Xx+c{t),t > 0} then conditionally to {S'^+'= = O}, M is almost 
surely finite. Thus there exists x > 1 such that 

P\M> x\S^+'^ = 0]< e'^iii^, 
L - ' J l-p(A + c) 



which implies 



P [5^+^ = 0|T^= < 00] = P [M > x\S'^^ = 0] ^if 3 



P [M > x] 



< e'. 



Furthermore for all e G [0, c], 

r , X4_ 1 P \S^^" = Ol - 1 



P = 0] - ^ 1 (16) 

- P[T^^<oo] + ' 



<e'. 
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Now let to > such that 



sup Xx{t) > x\T^ < oo 
te[o,io] 



We can find < c' < c such that 



> 1 - e'. 



(17) 



E (^o+^'(0 - V^{^) = 0, (Ar,VAA/t) n [0, to] = > 1 - e'. 

,«=i 

Conditionally to the event 

A := I sup > X, = 0, (AT^V^A/'a) H [0, to] 

we have almost surely T^~^^' < oo. Thus using ffTUj) . ffTTl) and ffTSj) 



(18) 



> P 



T^^"' <oo\A 



P [A\T^' < oo] = P [A|T,^ < oo] 



> (1 - e') > 1 - 2e'. 

Finally we get 

p{\) - p{\ + c') = P \S^^'' = 0, 5^ 
< p = o|T, 

<3e'. 

For the left continuity we do the same reasoning. 

Now we prove Theorem 14.31 in the case without interaction. 



A+c' 



< OO 



J- rv. 



A+c' 



oolT^ < oo 



□ 



Proposition 4.5. ^ and 



inf 



both have almost sure limits as t goes to infinity and 



1,. K 



■inf 



lim - lim 



A(l-A)p(A). 



e->-0 e t^oo t 

Proof. The convergences to almost sure limits are a consequence of Proposition 13.41 In- 
deed, for example Nt counts the number of elements of for which rit{l) 7^ 1 and 
r)t{2) = 1 minus the number of elements of JV^ for which rit{l) = 2. 
Then we have: 



inf 



t 



(19) 



i=l 



i=l 



Furthermore, since E [Si] = p{X) for alH > 1, the expectation of the right hand side of 
(fT9l) converges to A(l — A)p(A)e. But if we denote (tn)n>i the successive times for which 
the -th particle is exactly the n-th particle which will survive, then we have: 



at; 



inf 



-1 



i=l 
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Thus if := sup{n > I : tn < t} then, since ^ converges almost surely to (lim^. 
we have almost surely: 



inf 



-1 



lim 

t—*oo 



T\Tinf 

fit- 

t 



0. 



(20) 



On the other side, Ut is exactly the number at time t of particles which are already born 



and which will survive, i.e., rit = Yli=i almost surely. This shows that the right hand 



inf 



side of ffTI?|) converges also almost surely to A(l — A)p(A)e and fpUl) implies that 
converges to the same limit. □ 

4.2.2 Interaction implies survival 

The following lemma states that if a second class particle goes far enough, then it survives 
with high probability. 

Lemma 4.6. For all e' > 0, there exists xq (depending only on A and e') such that if c 
is small enough, then for all i > 1 

P [T^' <oo;3t> 0, X,{t) > xo] < e . 

Proof. We start by proving the same result for X**^-^. Let M := sup t > o|. 

Conditionally on < oo|, M is almost surely finite, thus we can choose Xq such 

that 



M > xom'''"^ < oo 



Hence 



7;^''"^<oo;3t>0,Xr-^(t) >a:o 



inf. 



M > xoIT:'''"^ < oo 



s,inf 



< OO 



< 



Furthermore, since the law of X*"-^ is the same for alH > 1 (because there is no interaction 
among supplementary particles), we can choose the same xq for alH > 1. Then we have 



P [T^' < oo; 3t > 0, Xi{t) > Xo] < P T^'"'"^ < oo; 3t > 0, -^(t) > xq 



Anf j 



+ P 



+ SI 



inf 



< e' 



if c is small enough. 



□ 



For X > 1, let Ta; := inf {t > : Xj(t) = x} (we omit the dependance on i in the 
notation because there will be no possible confusion). We can deduce from this lemma a 
stronger form of the same estimate: 
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Corollary 4.7. Let a; > 1. For all e' > 0, there exists xi depending only on X, e' and x 
such that if c is small enough then for all i > 1 

P [T,, <oo-3t> T,„X,{t) <x]< e'. 

Proof. We will use the same method as in Lemma [4. 1[ Let Et be the following event on 
the Poisson point processes of the Harris system during the time space [t,t + 1]: 

• one first class particle enters in site 1 and moves to site x; 

• then one first class particle enters and moves to site a; — 1; 

• we continue in the same way until x first class particles are entered in the system 
and they have moved until that the box {1, . . . , x} is full; 

• finally we impose that A4- fl [t, t + 1] = 0. 

Then ^^.(A) := P [Ef] depends only on A and x, is positive and under this event every 
second class particle which was in the box {1, . . . , x} at time t is died at time t + 1. 

Now let Xi be given by Lemma such that P [T/ < oo, 3t > 0, Xi{t) > xi] < e'g^(A) 
and define T+ := inf {t > T^^ : Xi{t) = x}. Then we have P [Tf < oo\T^ < c>o] > ^^.(A). 
This implies 

P [T,, < cx), 3t > T,,,X,{t) <x] = P \T+ < ool = ^^JJ ^ "^'.Jj" ^ < e'. 
L ^1 , _ XI, *v ; - J I X J p [7^^ < oo|T+ < oo] 

□ 

The next lemma states that if we fix x > 1, then the probability that two second class 
particles meet in the box {1, . . . , x} goes to with e. 

Lemma 4.8. Let Tj+i^, be the first time at which the {i + l)-th second class particle tries 
to jump on the site occupied by the i-th second class particle. Then for all fixed x > 1, 

P [Tj+i^j < oo, Xj(Tj+i^j) < x] — ^ 0, uniformly in i. 

Proof. Fix e' > and let xi and < cq < | — A be given by Corollary 14.71 such that 

P [T,., < oo, 3t > T^„Xi{t) <x]< e', for all e < Cq. (21) 
Then xi and cq depend only on A and e' (and x). We have: 



P[3s > t,Xi{s) e {l,...,x}] < P[3s > t,Xi{s) G {l,...,x},T,, < t] 

+ P[X,(t)>l,T,., >t], 
<P[T,., <oo,3s>T,„Xi{s) <x] 
+ P[Xi{s) G {l,...,Xi},VsG [0,t]]. 



(22) 



As in Lemma [4. 11 we have 
P[X,(s) G {l,...,Xi},VsG [0,t+l]] < (l-g,.,(A))P[X,(s) G {1, . . . , xj , Vs G [0,t]] 
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which imphes the existence of a constant C > depending only on A and e' such that 

P[X,(s)e{l,...,a;i},VsG[0,t]]<e-^*. 

Finally, using (!2T]) and (l22l) . there exists some deterministic to >0, depending only on A 
and e', such that 

P[3s>t,Xi{s) e {l,...,a;}] < 2e', 

for all t > to and e < cq. 

Besides, if we define as the time elapsed between Tf and the first jumping time 
of J\f2 greater than T/, then is an exponential random variable with parameter e 
independent of the trajectory of Xj. As a consequence, we have 

P [T,+i^, < oo, X,(T,+i^,) <x]<P[3t> T„ G {1, . . . , x}] , 

< P[3t > r„X,(t) G {!,..., x}, T, > to] 

+ P [Te < to] , 

Finally we have P [Tj+i^j < oo, Xj(Tj_|_i^j) < x] — > uniformly in i. □ 

Now we are able to prove that when a second class particle meets another one, both 
survive with a probability going to 1 as e goes to 0. 

Corollary 4.9. 

P [Tj+i^j < oo, T^,^ < ool — > 0, uniformly in i. (23) 

e— >0 

Proof. Fix e' > and let xq be given by lemma 14.61 We have 

P [Ti+i^i < oo,T-_^^ < oo] = P [Tj+i^i < 00,7^^+1 < oo, Xi(Ti+i_i) < xq] 

+ P [Ti+i^j < oo,7;^+i < oo,Xi+i(Ti+i^j) > xo] , 
< P [Tj+i^i < oo,Xi(Ti+i^j) < Xo] 
+ P [Tl^, < oo, 3t > 0, Xi+i(t) > Xo] , 

<26', 

if e is small enough. □ 

4.2.3 The proof of Theorem 14731 

Fix e' > and use fl23l) to find e > small enough to have 

P < oo,7;^+i < oo] <e'. 

We have already seen that both ^ and — ^ converge to almost sure limits and that 
- limt^oo converges almost surely to A(l — A)p(A) as e goes to 0. We also have 
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limj^oo ^ < A (1 — A) e, where is the number of second class particles entered in the 
system up to time t. Thus if we define 



r := inf <; t > : Vs > t, ^ < (A (1 - A) + 1) e 

s 



E 



rinf 



Ir 



<t 



(A(l-A)+l)et 



then r is almost surely finite and Nl"''^ — Nt = l^in/^^^^^ s {t)=o '^^ich implies 



< p sl-f{t) = l,S.{t) = 0,r<t 



=1 



< 



^ (A(l-A)+l)ei 



(24) 



i=2 

< (A(l-A) + l)ee', 



and, by dominated convergence theorem, the left hand side of fl24l) converges to lim^. 
lim^^oo ^ as t goes to infinity. Hence, dividing by e, we get 



inf 



1 N, 

< A (1 - A)p(A) - - lim — < (A (1 - A) + 1) e'. 



Since e' was arbitrary we can conclude. 
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